Revista Ingenierias Universidad de Medellin 5(9): 59-74 julio-diciembre de 2006 


Deficiencias, limitaciones, VENTAJAS y DESVENTAJAS 
DE LAS METODOLOGIAS DE ANALISIS SISMICO NO LINEAL 

Mauricio Alexander Mora*, Jesus Daniel Villalba** y Esperanza Maldonado*** 

Recibido: 27/07/2006 
Aceptado: 27/09/2006 


RESUMEN 

El presente trabajo esta dirigido a la comparacion de las metodologias de analisis sismico: dinamico 
cronologico no lineal, pushover tradicional, y dos procedimientos modales. La comparacion se realizo 
a partir del estudio de un edificio de mediana altura, aporticado, de concrete reforzado, representativo 
de las edificaciones tipicamente construidas en Bucaramanga. La estructura fue sometida a uno de los 
sismos definidos por el estudio de microzonificacion de la ciudad, escalado a diferentes valores de ace- 
leracion pico efectiva. Los resultados obtenidos fueron comparados en terminos de desplazamientos y 
derivas, con lo cual se pudieron detectar algunas deficiencias, limitaciones, ventajas y desventajas de las 
metodologias utilizadas, concluyendo que las metodologias modales son una buena aproximacion para 
la determinacion del comportamiento sismico de estructuras de concrete reforzado. 

Palabras clave: analisis sismico inelastico, analisis dinamico cronologico no lineal, analisis estatico no 
lineal, pushover tradicional, pushover modal. 

ABSTRACT 

This research presents the comparison of different non linear analysis methodologies: the non linear 
response history analysis, the traditional pushover methodology, and two modal pushover methodolo- 
gies. A medium rise reinforced concrete building, with moment resisting space frame, representative of 
the typical buildings in Bucaramanga, is examined using one of the earthquakes in the microzonation 
study of the city, which was scaled to different values of effective peak acceleration. The comparison of 
the results (floor displacements and interstory drifts) permitted to find deficiencies, limitations, advan- 
tages and disadvantages of the methodologies used. It was found that Modal Methodologies are a good 
approach to assess the seismic performance of reinforced concrete buildings. 
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INTRODUCCION 

El modelamiento del comportamiento de las 
estructuras ante eventos sismicos es complejo y 
la comunidad cientifica aun no ha definido un 
unico procedimiento para ello. 

En la actualidad este fenomeno puede ser modela- 
do de dos formas. La primera como un evento de 
naturaleza dinamica, para lo que se utiliza Analisis 
Dinamico Cronologico No Lineal (Non Linear Res- 
ponse History Analysis, NLRHA), el cual permite 
conocer la variacion de cualquier respuesta de la 
edificacion (desplazamientos de piso, derivas, fuerza 
cortante, etc.) en el tiempo. Asimismo, posee una 
base fisica rigurosa y es complicado de realizar. De 
otra parte, el efecto del sismo en las edificaciones 
puede modelarse como un evento estatico, median- 
te la utilizacion del procedimiento conocido como 
Analisis Estatico No Lineal o Pushover, descrito en 
algunos documentos como LEMA-273 o ATC-40. A 
pesar de tener algunas limitaciones desde el punto 
de vista metodologico y del alcance y exactitud de 
sus resultados, este procedimiento ha tenido una 
mayor aceptacion. 

El Analisis Pushover Tradicional (LEMA-273, ATC- 
40) solo tiene en cuenta el aporte del primer modo 
de vibracion, considerado imperante por poseer, 
generalmente, el mayor factor de participacion 
modal. Un aspecto importante en la aplicacion de 
esta metodologia es la construccion de la curva de 
capacidad de la edificacion. Esta curva se obtiene 
empujando el edificio horizontalmente mediante 
una distribucion de carga predefinida, la cual se in- 
crementa paso a paso hasta alcanzar un valor maximo 
de desplazamiento en el ultimo piso, o el colapso de 
la estructura. Los documentos guia no especifican 
un tipo de distribucion de carga a utilizar, sino que 
permiten cierta libertad. Lo anterior se convierte en 
una de las principales deficiencias de la metodologia, 
dado que la curva de capacidad del edificio es sensible 
a la distribucion de carga utilizada. 

Una posible solucion a las limitaciones de la tecnica 
pushover tradicional, antes mencionadas, es la apli- 


cacion de la metodologia pushover modal, en la cual 
se considera la participacion de varios modos. 

El principal objetivo de esta investigation es pre- 
sentar la metodologia de analisis modal, como una 
alternativa al riguroso analisis dinamico cronologi- 
co no lineal, conservando la facilidad de aplicacion 
de un analisis pushover tradicional y superando 
algunas de sus deficiencias y limitaciones. Lo ante- 
rior se lograra mediante la aplicacion de cada una 
de las anteriores metodologias y la comparacion 
de los resultados obtenidos en terminos de des- 
plazamientos y derivas de piso. 

DESCRIPCION 

DE LAS METODOLOGIAS 

Analisis dinamico cronologico no lineal 
(NLRHA, Nonlinear Response History 
Analysis) 

La esencia de la action sismica es la naturaleza di- 
namica de la carga, lo que conlleva a la variabilidad 
en el tiempo de las fuerzas internas del sistema y 
su estado (definido en terminos de desplazamien- 
tos y deformaciones). Por lo tanto, mientras un 
problema estatico tiene una unica solucion inde- 
pendiente del tiempo, la solucion de un problema 
dinamico requiere la descripcion del sistema en 
todos los instantes dentro del periodo en estudio, 
lo cual conlleva a grandes esfuerzos analiticos y 
numericos que pueden, a su vez, implicar el uso 
de una herramienta computacional. 

De otra parte, el analisis debe estar basado en una 
apropiada representation del movimiento del terreno 
y en principios aceptados por la dinamica estructu- 
ral. Para sistemas de multiples grados de libertad la 
ecuacion dinamica se puede expresar como: 

MG + CG-KU — MiU g (f) Ec. 1 

En donde M, C y K son las matrices de masa, amor 
tiguamiento y rigidez, respectivamente, U es el vector 
de desplazamientos, V es el vector de velocidad, U 
es el vector de aceleraciones, 0 (t) la historia de ace- 
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leraciones en el tiempo en la base del edificio, y cada 
elemento del vector de influencia X es igual a 1. 

El vector U esta relacionado con el numero de grados 
de libertad o posibilidades de movimiento de cada 
una de las masas concentradas. Generalmente, las 
masas se concentran en cada piso y se supone que no 
se pueden desplazar ni vertical ni rotacionalmente, 
permitiendo solo que se desplacen horizontalmente. 
Con estas suposiciones el sistema se reduce a dos 
posibilidades de movimiento o grados de libertad 
por masa: uno en el eje x y otro en el eje y. 

Si el anterior sistema de ecuaciones (Ec. 1) se plan- 
tea en una direccion, se tendra solo un grado de 
libertad para cada masa, y como estas se concentran 
en cada piso, el total de desplazamientos es igual a 
N, donde N es el numero de pisos; por consiguien- 
te, el vector U tiene tantos elementos como pisos 
tiene la edificacion. Igual sucede con las matrices 
M, C y K, las cuales tienen NxN elementos. 

La solucion de este sistema de ecuaciones es parti- 
cularmente complicada porque, si bien, la matriz 
de masa de la edificacion, M, es constante en el 
tiempo, las matrices de amortiguamiento y rigidez, 
C y K, respectivamente, cambian cuando se supera 
el limite de comportamiento elastico. 

Sin embargo, reconstruir las matrices de rigidez y 
amortiguamiento cada vez que un elemento (viga, 
columna, muro) supera la fluencia no es el unico 
problema. La parte derecha de la ecuacion (1) re- 
presenta las fuerzas externas a las que se somete la 
edificacion, siendo UJt) la aceleracion en la base 
del edificio en el tiempo, es decir, el acelerograma, 
de manera que la ecuacion de movimiento es dife- 
rente para cada valor de aceleracion. 

En el presente estudio, el Analisis Dinamico Crono- 
logico No Lineal (NLRHA) se lleva a cabo utilizando 
el software IDARC 2D version 6 (Inelastic Damage 
Analysis of Reinforced Concrete Structures; State 
University of New York at Buffalo, Department of 
Civil Engineering, Julio 2004) y se asume que los 
resultados obtenidos mediante este tipo de analisis 
corresponden al punto de referencia, es decir, que 


estos serviran para compararlos con los obtenidos 
por los metodos de analisis no lineal. 

Analisis estatico no lineal-pushover tradicional 

El Analisis estatico no lineal o analisis pushover 
es otro tipo de metodologia para determinar la 
respuesta no lineal de la estructura. A diferencia 
del analisis dinamico cronologico no lineal, en 
la metodologia pushover se calculan los valores 
maximos de cualquier respuesta de la edificacion, 
y no su historia en el tiempo. 

El analisis tiene como base dos conceptos pri- 
mordiales: la capacidad y la demanda. La primera 
puede entenderse como una caracteristica propia 
de la edificacion que depende de factores como la 
geometria de los elementos, la cantidad de refuer- 
zo, las propiedades de los materiales como rigidez 
y ductilidad, entre otras. La demanda depende de 
un sismo en particular, representado mediante 
un acelerograma o un espectro de diseno, y hace 
referencia a las fuerzas y deformaciones impuestas 
por este. De manera que la demanda, a diferencia 
de la capacidad, no es en teoria un valor constante, 
pues depende del conjunto de fuerzas externas o 
aceleraciones a las que se somete la edificacion. 

Para evaluar la capacidad de la edificacion se debe 
construir la curva de capacidad, la cual se obtiene 
graficando los valores maximos de cortante basal y 
desplazamiento en el ultimo piso. El procedimiento 
para construir la curva de capacidad consiste en 
aplicar cargas laterales monotonicas, las cuales se van 
incrementando hasta llevar la estructura al colapso. 
En la figura 1 se muestra un esquema de las fuerzas 
que se aplican por piso a la estructura; la sumatoria 
de todas las fuerzas F. es igual al cortante basal V. 

La edificacion responde a la aplicacion de las cargas, 
desplazandose lateralmente hasta que alcanza un maxi- 
mo u_. El par de valores uyV representa un punto 
en la curva de capacidad de la edificacion, de manera 
que para determinar otros puntos se carga la estructura 
nuevamente con fuerzas mayores y, asi, hasta generar 
un mecanismo de colapso en la estructura. 
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Figura 1, Procedimiento pushover para determinar la curva de capacidad de la edificacion 


La demanda varla dependiendo del sismo al que 
se somete la edificacion y se evalua mediante 
espectros de respuesta (o diseno), El concepto 
de espectro de respuesta aparece como un medio 
conveniente para resumir la respuesta maxima de 
sistemas lineales y no lineales de un grado de liber- 
tad (SDF) a un movimiento sismico en particular, 
Igualmente, permite aplicar el conocimiento de 
la dinamica estructural al diseno de estructuras y 
al desarrollo de los requerimientos de resistencia 
a fuerzas sismicas consignados en los codigos 
(Chopra, 2001a), 

En la literatura pueden encontrarse articulos de 
varios expertos, relacionados con este tipo de ana- 
lisis, tales como Krawinkler y Seneviratna (1998), 
Mwafy y Elnashai (2001), Inel et al (2003), Lopez 
(2004), Sierra (2004), entre otros, 

Analisis no lineal modal 

En el ano 2001, Anil K, Chopra y Rakesh K, Goel 
publicaron el trabajo titulado “A modal pusho- 
ver procedure to estimate seismic demands for 
buildings: Theory and preliminary evaluation”, 
el cual es una combinacion de las metodologias 
presentadas anteriormente (Analisis dinamico 
cronologico no lineal y estatico no lineal), Ini- 
cialmente, se implementa un Analisis estatico no 
lineal (pushover) para evaluar la capacidad por 


modo; posteriormente, todo el sistema de N grados 
de libertad (MDF) se reduce a N sistemas de un 
solo grado de libertad (SDF), que se someten a la 
accion sismica, para, finalmente, combinar todos 
los aportes modales y encontrar una aproximacion 
a la respuesta total, 

El procedimiento tambien se presenta como una 
solucion a las deficiencias encontradas en el Ana- 
lisis pushover tradicional, ya que se define una 
sola distribution de carga horizontal y toma en 
cuenta el aporte de todos los modos de vibration 
a la respuesta total de la estructura, 

Inicialmente, para aplicar este metodo, es nece- 
sario desarrollar un procedimiento de Analisis 
modal desacoplado de la historia de respuesta 
(UMRHA), que es una buena aproximacion al 
Analisis dinamico cronologico no lineal (NLRHA) 
o Analisis dinamico, para sistemas de multiples 
grados de libertad (MDOF), 

Analisis modal desacoplado de la historia 
de respuesta (UMRHA, Uncoupled Modal 
Response History Analysis) 

La respuesta inelastica de un edificio de N-pisos 
con simetria en planta en los dos ejes ortogonales x 
y y, cuando se somete a un movimiento sismico en 
la direction de uno de sus ejes de simetria (x o y), 
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se puede determinar como funcion del tiempo me- 
diante el analisis UMRHA (Chopra y Goel, 2001). 
El procedimiento a seguir parte de determinar la 
frecuencia natural <d y los modos de vibracion (b 
bajo comportamiento lineal; para el i-esimo modo 
se elabora la curva de capacidad utilizando una sola 
forma de distribution de carga s n \ que depende de 
la matriz de masa y de la forma modal. La curva 
obtenida es idealizada como una relacion bilineal y 
luego convertida a la relacion F sr /L n - D n , utilizando 
las siguientes ecuaciones: 


ij 

>rjL' 
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Ec. 2 


Siendo V bny el cortante en el punto de fluencia; 
M „ la masa modal efectiva; D es el desplaza- 
miento de fluencia del sistema de un grado de 
libertad equivalente del MDF para el modo n; 
u rny el desplazamiento de fluencia en el ultimo 
piso; F /L es la fuerza resistente; el factor de 
participation modal. 

A su vez, se calcula la historia de deformacion, DJt), 
y la seudo-aceleracion, A n (t), para el i-esimo modo 
de vibracion del sistema inelastico de un grado 
de libertad asumiendo que la relacion fuerza-de- 
formation es la relacion F /L - D . Con ello, se 

sn n n 

determina la historia de varios tipos de respuesta 
como derivas de piso, fuerza cortante en la base, 
etcetera, mediante las ecuaciones: 

M 0 = A <*>; rjt) = rfA„ U) Ec. 3 

En donde ujt) representa la historia en el tiempo de 
los desplazamientos de piso y rjt) es la variacion en 
el tiempo de cualquier tipo de respuesta, teniendo 
en cuenta que r st n corresponde a la respuesta modal 
estatica, es decir, el valor estatico de r (respuesta) 
debido a las fuerzas externas s* . 

n 

El procedimiento anterior se realiza para tantos 
modos como sean necesario para obtener una 
buena exactitud. Por lo general, los primeros dos 
o tres modos son suficientes. 


Finalmente, se combinan las respuestas modales 
para determinar la total utilizando las siguientes 
ecuaciones: 


N N 

uit) = X ^) = X T r$n D n (0 


n=\ n=\ 


.V j¥ 




rj=l hj=I 


Ec. 


4 


A partir del resultado anterior se puede determinar 
el valor pico, r y de la respuesta total 

Analisis pushover modal 
(MPA, Modal Pushover Analysis) 

Chopra y Goel (2001), propusieron el Analisis 
pushover modal (MPA) como una herramienta 
para determinar la respuesta inelastica pico de un 
edificio sometido a un sismo. 

En este metodo se desarrollan los primeros pasos 
del analisis modal desacoplado (UMRHA) hasta 
convertir la curva de capacidad bilineal en la 
relacion F /L - D . Luego se determina la defer- 
macion pico, D , del i-esimo modo de vibracion del 
sistema inelastico de un grado de libertad (SDF) 
utilizando la relacion fuerza-deformacion. Luego, 
se calcula el desplazamiento maximo en el ultimo 
piso u asociado al vesimo modo de vibracion del 
sistema inelastico SDF utilizando la ecuacion: 

u rno = r A D n Ec - 5 

Cuando se alcance la deformacion u , se extrae 

r no 

de la base de datos pushover los valores de otras 
respuestas, r no . Se repite lo anterior para tantos 
modos como sean necesarios, hasta obtener sufi- 
ciente exactitud. Al igual que para la metodologia 
UMRHA la utilizacion de los dos o tres primeros 
modos es suficiente, por lo general. 
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Finalmente, se determina la respuesta total utili- 
zando la regia de combinacion SRSS (raiz cuadrada 
de la suma de cuadrados): 


o 



Ec, 6 


APLICACION 

DE LAS METODOLOGIAS 


Seleccion de la estructura 

La estructura utilizada corresponde a un edificio 
de concreto reforzado, de 7 pisos, regular en 
planta y altura, representativo de las edificaciones 
construidas en el municipio de Bucaramanga, El 
sistema de resistencia sismica esta conformado 


por porticos ortogonales entre si, Cuenta con una 
planta de 1L6 m * 10,5 m, y una altura de 17,7 m, 
En la figura 2 se presenta una vista en planta de los 
niveles L6 y la elevacion de los porticos en cada 
direccion, Los porticos son iguales en cada direc- 
cion: 3 luces de 4.4, 2,45 y 4.75 m en la direccion 
este-oeste (E-W: direccion x) y 2 luces y 2 voladizos 
exteriores de 4.15 y 1,1 m respectivamente en la 
direccion norte-sur (N-S: direccion y), 

El sistema de entrepiso esta conformado por placas 
rigidas aligeradas apoyadas en la direccion esteoeste 
(E-W), La planta del nivel 7 es igual a la de los demas 
niveles, excepto por los voladizos. Las columnas son 
de seccion cuadrada, de 0,5 X 0,5 m y las vigas rectam 
gulares de 0,35 x 0,45 m. La resistencia del concreto 
utilizado en placas, vigas y columnas es de 21 MPa, 
y la del acero de refuerzo tanto longitudinal como 
transversal es de 420 MPa (Mora y Villalba, 2005), 



r 






Figura 2* (a) Vista en planta de los pisos L6; (b) Porticos A, B, C; (c) Porticos 1, 2, 3 y 4; y Refuerzo 
longitudinal y transversal: (d) columnas y (e) vigas, Todas las medidas en metros 
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Definicion de la accion sismica 

Para la aplicacion de las metodologias se utilizo el 
sismo de Coalinga, el cual garantiza que la estruc- 
tura se comporte inelasticamente para el intervalo 
de aceleraciones que se va a trabajar. El registro 
fue tornado en la estacion USGS 1608 Oil Fields 
Fire Station a las 02:39 del 22 de julio de 1983, 
presento una magnitud de 5.77, profundidad de 7.4 
km y aceleracion pico efectiva de 0.213g. Este sismo 
es uno de los definido por INGEOMINAS para 
la realization del “Estudio de microzonificacion 
sismica del area metropolitana de Bucaramanga” 
(Ingeominas, 2001). En la figura 3 se presenta el 
registro de los 16 primeros segundos del sismo, 
dentro de los cuales se presento la aceleracion 
pico efectiva. 



Figura 3* Sismo Coalinga 


Analisis estatico no lineal- pushover 
tradicional 

Con el fin de observar la influencia de la distri- 
bution lateral de fuerzas en la construccion de la 
curva de capacidad, esta se calculo para la direction 
x del edificio en estudio, tomando una distribucion 
uniforme, una distribucion ELF (equivalent lateral 
force) y una triangular invertida, con los valores 
que se muestran en la figura 4* 


77*i m — 

949KN — 

949KN — 

949 KN — 

1-452 KN 

iSOO KN 

1 226 KW 

95? m 

&94HM — 
W HW - 
1 99 KN 



Figura 4 ♦ Distribucion lateral de fuerzas en 
la construccion de la curva de capacidad. (a) 
Distribucion uniforme; (b) Distribucion ELF; (c) 
Distribucion triangular invertida o lineal 

Para la construccion de las curvas de capacidad 
se utilizo el software IDARC 2D; los resultados 
se presentan en la figura 5. Estos muestran que 
existe una considerable variacion en funcion de 
la distribucion de carga horizontal utilizada. Las 
distribuciones ELF y triangular invertida son muy 
similares, pero difieren de la distribucion uniforme, 
incluso en la rigidez inicial. 
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DESPLAZAMIENTO ULTIMO PISO [cm] 

Figura 5* Curvas de capacidad utilizando 
diferentes distribuciones de carga horizontal 

El siguiente paso consiste en determinar la demanda 
sismica, En primer lugar, es necesario definir nume- 
ricamente el movimiento sismico, primero como 
acelerogramay despues como espectro de demanda, 
A manera de calculo tipo se utilizo el sismo Coalinga, 
escalado uniformemente a 0,5g (figura 3), 

El punto maximo hasta el que es llevada la edificacion 
se conoce como punto de comportamiento, y para en- 
contrarlo se han propuesto diferentes metodologias, 
En esta investigation se utilizo el procedimiento A del 
ATC-40, que es basicamente un proceso de prueba y 
error: una vez se tienen en una misma grafica (forma- 
to ADRS) el espectro de capacidad de la edificacion 
y el espectro elastico de demanda (amortiguamiento 
= 5 %), se supone un punto de comportamiento 
cualquiera, luego se construye una representation 
bilineal de la capacidad, que debe tener la misma 
pendiente inicial que el espectro, y la misma area 
bajo la curva hasta el punto que se asume es el punto 
de comportamiento, Posteriormente, se determina 
el amortiguamiento desarrollado por la estructura y 
se construye el espectro de demanda para este valor 
de amortiguamiento. Si el espectro de demanda 
calculado intercepta al espectro de capacidad en el 
punto previamente definido, entonces la suposicion 
fue correcta y ese es el punto de comportamiento. 
Si los puntos no coinciden, es necesario suponer un 
nuevo punto de comportamiento y repetir el proceso 
tantas veces como sea necesario hasta que los dos 
espectros se intercepten en ese punto, En la figura 6 


se presenta el punto de comportamiento para cada 
distribucion de carga y el espectro elastico para el 
sismo de Coalinga, 


Una vez conocido el punto de comportamiento es 
posible determinar casi cualquier valor de respuesta 
maxima de la edificacion, En la figura 7(a) puede ver- 
se el perfil de desplazamientos de piso maximos como 
porcentaje de la altura de la edificacion (17,7 m), 
obtenido mediante cada una de las tres distribuciones 
horizontales de carga y el perfil asumido Verdadero” 
mediante la aplicacion del Analisis dinamico crono- 
logico no lineal, De igual forma, en la figura 7(b) se 
presenta el perfil de maximas derivas de entrepiso 
como porcentaje de la altura de cada piso. 



Figura 6 ♦ Punto de comportamiento: (a) 
Distribucion uniforme (Peq = 28,74); 

(b) Distribucion ELF (Peq = 30,51) y (c) 
Distribucion triangular invertida (Peq = 30,49) 
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Figura 7* Perfiles de (a) desplazamiento de piso maximos y (b) maximas derivas de entrepiso, 
normalizados por la altura total de la edificacion 


Con el fin de observar la influencia de la intensidad 
del movimiento de tierra, el sismo Coalinga se esca- 
1 6 para que alcanzara una aceleracion pico efectiva 
igual a 0.2g, 0.3g, 0.4g, 0.6g y 0.7g. Para cada valor 
se construyo el espectro de respuesta y se aplico el 
mismo procedimiento descrito anteriormente. 

En las figuras 8 y 9 se presentan resumenes con 
los porcentajes de error en los desplazamientos de 


piso maximos y en los valores pico de las derivas de 
entrepiso, obtenidos para la distribucion uniforme, 
ELF y triangular, respectivamente. Se asume que 
los resultados obtenidos por el NLHRA correspond 
den a los mas cercanos a la realidad, por lo cual, la 
validacion y comparacion de los resultados se hizo 
en funcion de ellos. 




4* Jib (M\ liifTU 



Figura 8* Error en los maximos desplazamientos de piso: (a) Distribucion uniforme; (b) Distribucion 

ELF; (c) Distribucion triangular 
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Figura 9* Error en las derivas de entrepiso maximas: (a) Distribucion uniforme; (b) Distribucion ELF; 

(c) Distribucion triangular 


En el caso de los desplazamientos, puede observar- 
se que no siempre la misma distribucion horizontal 
de carga genera los mayores o menores porcentajes 
de error. Por ejemplo, para 0.5g, los porcentajes de 
error mas grandes se obtienen para la distribucion 
uniforme, pero para 0.3g se obtienen los menores. 
Con las derivas pasa lo contrario los porcentajes de 
error mas grandes, sin importar el factor de escala, 
se presentan para la distribucion uniforme. Las 
distribuciones ELF y triangular generan errores 
muy similares, debido a la similitud de la curva 
de capacidad obtenida, ver Figura 5. 

Los errores, en algunos casos mayores que el 50%, 
se pueden relacionar con: 

1 . Las suposiciones inherentes a la teoria utilizada 
para obtener la expresion que define el calculo 
del amortiguamiento equivalente (p ). 

2. La idealizacion del espectro de capacidad 
como una relacion bilineal. 

3. El espectro de respuesta se obtiene para 
un sistema de un grado de libertad (SDF), 
y luego se aplica a un sistema de multiples 
grados de libertad (MDF). 

4* La suma de errores, primero al igualar las 
areas bajo el espectro de capacidad real y la 


representacion bilineal y, luego, al igualar el 
punto de comportamiento supuesto con el 
punto de corte de los espectros de capacidad 
(real) y demanda. 

5. El procedimiento solo tiene en cuenta el 
aporte del primer modo de vibracion. 

6. La distribucion de carga horizontal seleccio- 
nada. 

Analisis estatico no lineal modal 

La respuesta del edificio de siete pisos se deter- 
mino por medio de los dos metodos modales 
inelasticos aproximados UMRHA (Analisis 
modal desacoplado de la historia de respuesta) 
y MPA (Analisis pushover modal) y, se comparo 
con los resultados del riguroso procedimiento 
NLRHA utilizando el software IDARC-2D. Para 
garantizar que la edificacion respondiera dentro 
del rango inelastico, inicialmente se utilizo el 
acelerograma del sismo Coalinga escalado a 0.6g. 
En la figura 10 se presentan las propiedades di- 
namicas de los tres primeros modos de vibracion 
del edificio en la direccion x y la distribucion 
lateral de fuerzas: 
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Figura 10* (a) Forma y periodo de los tres primeros modos 
y (b) Distribucion de fuerzas horizontales, en KN 


Las curvas de capacidad, para cada uno de los tres 
primeros modos de vibracion obtenidas utilizando 
la distribucion lateral de fuerzas de la figura 10, se 
presentan en la figura 11. 

Como se menciono anteriormente, el analisis 
UMRHA permite determinar una aproximacion a 
la historia en el tiempo de los desplazamientos de 



piso, o cualquier otra respuesta de la edificacion. 
En la figura 12 se presenta la variacion del despla- 
zamiento en el primer y ultimo piso, obtenida me- 
diante el analisis UMRHA y mediante el NLRHA. 
Puede verse que el analisis UMRHA constituye una 
buena aproximacion y, en algunos casos, puede 
utilizarse como alternativa, pues es mas facil de 
implementar que el analisis NLRHA. 
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Figura 11. Curvas de capacidad para los tres primeros modos de vibracion (direccion x) 
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Figura 12* Metodo UMRHA vs NLRHA. (a) Historia de desplazamientos del piso 7 y (b) Historia de 
desplazamientos del piso L 


En la figura 13(a) se muestra el perfil de desplazamientos de piso maximos, obtenidos mediante UMR- 
HA y NLRHA. El porcentaje de error mas grande ocurre en el piso 5, en donde se alcanza un 14*5% y 
el menor en el piso 2, con 3.27%* 




Figura 13* Metodo UMRHA vs NLRHA. (a) Perfil de desplazamientos maximos (Sismo Coalinga 
escalado a 0.6g) y (b) Error en los desplazamientos de piso para diferentes factores de escala 
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Con el fin de observar el cambio de los porcentajes 
de error en la medida en que se varia la magnitud del 
sismo, el mismo procedimiento se llevo a cabo para 
el sismo Coalinga escalado a O.lg, 0.2g, 0.3g, 0.4g, 
0.5g y 0.7g. Los resultados se presentan para cada 
piso en la figura 13(b). Como era de esperarse, los 
errores mas pequenos se presentan para los menores 
factores de escala (O.lgy 0.2g) pues la estructura aun 
no ha incursionado en el rango inelastico. Para fac- 
tores de escala intermedios (0.3g, 0.4gy 0.5g) se pre- 
sentan los porcentajes de error mas grandes debido 
a que corresponden a la zona de menor correlation 
entre la curva de capacidad real y la idealizada, que 
es la que en realidad se utiliza. El error tiende a 
disminuir cuando la estructura ha incursionado en 
gran medida en la fluencia (0.6g y 0.7g). 


Finalmente, los maximos desplazamientos de piso y 
maximas derivas de entrepiso del edificio aplicando 
el sismo Coalinga escalado a 0.6g en la direction x, 
se evaluaron mediante el Analisis pushover modal 
(MPA) tomando como base las curvas de capacidad 
por modo de la figura 11. En las figuras 14 y 15 
se presenta el perfil de desplazamientos y derivas 
maximas para cada piso, teniendo en cuenta uno, 
dos y tres modos, y el perfil calculado mediante 
NLRHA. Tambien se presenta el error teniendo en 
cuenta el primer modo, los dos primeros modos y, 
finalmente, los tres primeros modos de vibracion 
del edificio. El analisis de las graficas nos permite 
concluir que el primer modo de vibracion no es 
suficiente cuando se quiere implementar un ana- 
lisis modal. 



Figura 14* Metodo MPA vs. NLRHA. (a) Perfil de desplazamientos maximos; 
(b) % de error considerando uno, dos y tres modos 
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Figura 15. Metodo MPA vs. NLRHA. (a) Perfil de derivas maximas; 
(b) % de error considerando uno, dos y tres modos 
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El analisis MPA se implemento adicionalmente 
para el sismo Coalinga escalado a 0,lg, 0.2g, 
0,3g, 0,4g, 0,5g y 0,7g- Nuevamente se tuvo en 
cuenta el aporte de los tres primeros modos de 
vibracion y se determinaron los valores pico de 


desplazamientos de piso y derivas de entrepiso, 
Los errores relativos al analisis NLRHA se pre- 
sentan en la figura 16, encontrandose conclu- 
siones similares a las halladas para el metodo 
UMRHA, 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 


Factor de escala del sismo 


Factor de escala del sismo 


Figura 16 * Analisis Pushover Modal (MPA), Error como funcion del factor de escala del sismo; (a) 

Desplazamientos y (b) Derivas de entrepiso 


CONCLUSIONES 

El Analisis dinamico no lineal (NLRHA) es la 
metodologia de analisis sismico mas completa que 
existe en la actualidad, pues permite conocer la 
variacion en el tiempo de cualquier respuesta de 
la edificacion, Sin embargo, existen metodologias 
alternativas mas utilizadas, debido principalmente 
a que son mas sencillas de implementar, aunque 
menos exactas, Para determinar el grado de 
exactitud de estas, los desplazamientos de piso y 
derivas de entrepiso de un edificio de concrete 
reforzado, aporticado, de siete pisos, regular en 
planta y altura, sometido al sismo Coalinga, fueron 
comparados contra aquellos obtenidos mediante 
el Analisis dinamico no lineal, encontrandose que 
los mayores porcentajes de error se presentaron 
para el procedimiento pushover tradicionah 

Se encontro que la construccion de la curva de 
capacidad, necesaria para implementar un ana- 
lisis pushover, es muy sensible a la distribucion 
horizontal de carga, lo que en teoria no deberia 
suceder, ya que la capacidad de la edificacion es 
unica, Para las distribuciones ELF (equivalent late- 


ral force) y triangular las curvas son practicamente 
iguales, pero existe gran diferencia, incluso en la 
rigidez antes de la fluencia, con la curva obtenida 
mediante la distribucion uniforme, 

Tambien se puede concluir que la distribucion 
de carga que genera los errores mas grandes es la 
distribucion uniforme, Los errores disminuyen 
utilizando las otras dos distribuciones, sin embar- 
go, siguen siendo mayores que aquellos obtenidos 
mediante los procedimientos UMRHA y MPA, 
Las explicaciones son basicamente dos; primero, 
el procedimiento pushover tradicional solo tiene 
en cuenta el primer modo de vibracion, y segun- 
do, la incertidumbre en torno a la seleccion de 
la distribucion horizontal de carga con la que se 
construye la curva de capacidad, Puede verse que 
este factor tiene gran influencia en los resultados 
obtenidos, y no existe un criterio unificado para 
seleccionar solamente una, 

Los dos problemas mencionados en el parrafo an- 
terior se pueden solucionar con los procedimien- 
tos UMRHA y MPA, En los dos procedimientos 
modales se utiliza solamente una distribucion 


Revista Ingenierfas Universidad de Medellin 


Deficiencies, limitaciones, ventajas y desventajas ... 


horizontal de carga, proporcional a la masa y forma 
del modo, para construir la curva de capacidad; de 
esta manera se obtiene solo una curva para cada 
modo y se elimina la incertidumbre sobre cual dis- 
tribucion utilizar, presentada en el procedimiento 
pushover tradicional. 

Los errores obtenidos en las dos metodologias mo- 
dales (UMRHA y MPA) aumentan cuando la edifi- 
cacion ingresa en el rango inelastico, reflejando de 
esta forma la dificultad e incertidumbres existentes 
para tratar de predecir el comportamiento no lineal 
de las estructuras. Las causas de dichos errores son, 
en primer lugar, la falta de una teoria rigurosa, 
principalmente en el analisis UMRHA. Los dos 
procedimientos se basan en algunas suposiciones 
y extensiones del analisis elastico de edificaciones 
que no son del todo ciertas cuando el analisis es 
inelastico. En segundo lugar, se induce al error al 
idealizar la curva de capacidad como una relacion 
bilineal. Los resultados permitieron determinar 
que los mayores errores se cometen en la parte de 
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